arXiv: 1506.00805vl [physics.atom-ph] 2Jun2015 


Can one turn off Coulomb focusing? 

S.A. Berman 

School of Physics, Georgia Institute of Technology, 

Atlanta, Georgia 30332-0430, USA and 
Gentre de Physique Theorique, GNRS-Aix-Marseille Universite, 

Gampus de Luminy, 13009 Marseille, Franee 

C. Chandre 

Gentre de Physique Theorique, GNRS-Aix-Marseille Universite, 

Gampus de Luminy, 13009 Marseille, Franee 

T. Uzer 

School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332-0430, USA 

(Dated: June 3, 2015) 

Abstract 

We find that Coulomb focusing persists even when the Coulomb field is barely noticeable com¬ 
pared with the laser field. Delayed recollisions proliferate in this regime and bring back energy 
slightly above the 3.17 Up high-harmonic cutoff, in stark contradiction with the Strong Field Ap¬ 
proximation. We investigate the nonlinear-dynamical phase space structures which underlie this 
dynamics. It is found that the energetic delayed recollisions are organized by a reduced number of 
periodic orbits and their invariant manifolds. 
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I. INTRODUCTION 


The theoretical framework for electron dynamics processes in intense linearly polarized 
laser fields was established about twenty years ago and remains the state-of-the-art in strong 
field physics [1-3]. It centers on the “recollision” model which follows a three-step scenario: 
Electrons are first ionized, absorb energy while following the laser, and then are propelled 
back to the ionic core after the laser reverses direction, about half a laser cycle after ion¬ 
ization. The kinetic energy transfer in the core region causes ionization of more electrons 
(nonsequential multiple ionizations - NSMI) [4] or generates very high harmonics of the 
driving laser by high harmonic generation (HHG) [5]. Despite the wide-spread use of the 
three-step model, it is not widely appreciated that it usually takes more than one recollision 
for the ionized electron to transfer its energy to set the processes of multiple ionization or 
high-harmonic generation into motion. It is reasonable to expect that the spreading of the 
electronic wave packet would make these delayed collisions much less effective. Surprisingly, 
the opposite is the case, and it is due to a process called Coulomb focusing [6]. The ef¬ 
fectiveness of Coulomb focusing was demonstrated convincingly in the pioneering work of 
Brabec et al. [6], who showed that the nonsequential double ionization rate is enhanced by 
an order of magnitude by Coulomb focusing. 

Research in the last two decades has confirmed that Coulomb focusing (and the Coulomb 
field in general) is a key player in the recollision process [7-11]. But it seems that the 
opposite question -namely, when the Coulomb focusing does not matter- has not received 
the same attention. Perhaps the reason is that the answer seems obvious: the Coulomb 
field can be ignored when the laser field is strong. After all, this is one of the tenets of the 
Strong Field Approximation (SEA) [12] which ignores the Coulomb field when the electron 
is moving in the laser field. But how strong is “strong” for the laser field, and how weak 
does the Coulomb field have to be for us to neglect it altogether -i.e., is it possible to turn 
it off completely, as the three-step model would have us do? 

In fact, the Coulomb field and its focusing effect persist well beyond what would be 
expected based on the comparison of field strengths, and well beyond field strengths where 
there is no potential barrier to trigger tunneling ionization at peak field. In this range of 
intensity, contrary to SEA predictions, delayed recollisions continue to manifest themselves 
and bring back energy slightly above the 3.17 Up high-harmonic cutoff. The kinetic energy 
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brought back by delayed recollisions is compatible with the value provided in the SFA 
for immediate recollision (occurring half a laser cycle after ionization), but those energetic 
delayed recollisions do not exist in the SFA. 

Our paper is structured as follows: In Sec. II we begin by reviewing how the return energy 
is maximized in the SFA, and characterize the trajectories bringing maximum energy to the 
core. In Sec. Ill we investigate how the energy-maximizing classical recollision trajectories 
change when the Coulomb field is included. We show that delayed energy-maximizing tra¬ 
jectories, which are not present in the SFA, emerge. We provide a characterization of the 
delayed recollisions in phase space, and show that they are related to the invariant manifolds 
of particular periodic orbits. 


II. MAXIMIZING THE RETURN ENERGY IN THE STRONG FIELD APPROX¬ 
IMATION 

We begin with the derivation of some results on recollisions in the Strong Field Approx¬ 
imation (SFA). These results will form the basis for understanding the distinct effect of the 
Coulomb field on recollision trajectories. The SFA may be described by the one-dimensional 
classical Hamiltonian (in atomic units) 

H{x,p,t) = ^(^p+^cos{ujt + (l))^ , (1) 

where p is the canonical momentum, Eq is the electric field amplitude, uj is the laser fre¬ 
quency, and the dipole approximation for the laser field is used. In the context of the 
three-step model, this Hamiltonian describes the motion of the electron after its ionization, 
which we define as occurring at t = 0 without loss of generality (up to a change of phase 
(p). The three-step model does not specify the exact phase of the laser at the instant of 
ionization, so we allow this phase 0 to be a free parameter of the model. 

We write Hamiltonian (1) in the velocity gauge as opposed to the more common length 
gauge to make the translational invariance of the SFA apparent (i.e., independence of Hamil¬ 
tonian (1) with respect to x). Thus, it is clear that the momentum p is conserved, and as a 
consequence, the resulting system is integrable. It is convenient to rescale the variables so 
that X is in units of the quiver radius ^ p is in units of Eq/lo, and t is in radians, so 


3 



( 2 ) 


Hamiltonian (1) now looks like 

'H{x,p,t) = ]^(p + cos{t + (j))^ • 
We also note that this Hamiltonian has the symmetry 


(x,p,0) ^ (-X,-p,7r + (/)), 

which allows us to narrow our focus to initial conditions with 0 G [0, tt). 
From Eq. (2) we write down the solution for x(t), 

x(t) = Xq — sin (p + Pot + sin(t + 0 ), 


( 3 ) 


( 4 ) 


and p{t) = pq. Just as the three-step model does not exactly specify the laser phase p at 
ionization, it also does not specify the position of the electron Xq at ionization. Thus, we 
allow Xq to be a free parameter as well. Also, we note that Eq. (4) shows that Pq is the 
constant drift momentum of the electron. 

Equation (4) is used to compute the maximum return energy of the electron. First, it is 
assumed that the electron reaches the atomic core at time A, he. x{tr) = 0. By substituting 
this into Eq. (4), one obtains 


Po(trt P') 


sin p — sin (A + p) — Xq 


( 5 ) 


Now, one may fix Xq and maximize the energy at return, 


(j)) 


Po + cos(t^ 


2 

•> 


( 6 ) 


with respect to tr and p. Setting the derivatives equal to zero gives 


= (po + cos(A + P)J (dt^Po - sin(A + p)^ = 0, (7) 

d^e = + cos(A + P)^ [d^yPo - sin(A + P)^ = 0. (8) 

We discard the solution po + cos (A + ^) = 0 which corresponds to the zero-energy minima 

of Eq. (6). Thus, we get dt^po = d^po, and carrying out these derivatives using Eq. (5) 

immediately gives 

x{t = 0) = po J- cos 0 = 0, (9) 


which corresponds to zero initial momentum in the length gauge. Though this result has 
been known for xq = 0 , we have shown that it is independent of this condition and is valid 
for any Xq 7 ^ 0 . 
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In addition, Eq. (9) tells us that when the return energy is maximized, the initial energy 
is zero and thus minimized, implying that the change in the energy 

Ae = ^ (po + cos{tr + - (^po + cos 

is at a maximum as well. By rewriting Ae as 

Ae = (^Pq + - COs{tr + ((>) + - cos (f)^ ^ COs{tr + (f>) — cos (f)^ , (10) 

substituting Eq. (5) into Eq. (10), and using some trigonometric identities, it is shown that 


Ae{tr, 0) = f(tr) sin(t^ + 20) + g(tr, 0)a:o, 


where /(t^) = 
and g{tr, 0) = 


sin y ( 2 sin y — tr cos y 


2 sin sin(A + 


We first treat the case of the electron starting very close to the core, i.e. xq 
case Eq. (11) simplifies to 


( 11 ) 
( 12 ) 
(13) 
0. In that 


ASoitr, 0) = f{tr) Sm{tr + 20). (14) 

To compute the maximum of this function, we first note that all the 0 dependence is con¬ 
tained in the multiplicative sine term. Also, one may find a 0 such that sin(tr + 20) = ±1 
for any tr- Therefore, the absolute maximum of Asq occurs where |/(tr)| is at its absolute 
maximum for tr > 0, and sm{tr + 20) is the sign of f{tr) at the extremum. Then we obtain 
the absolute extremum of f(tr) in Eq. (12) by setting dt^f = 0 and using trigonometric 
double angle identities, yielding 

2 — 2tr sin tr + {tl — 2) cos tr = 0, (15) 

in agreement with the literature [2, 5, 13]. This equation’s first positive root, t*Q ~ 4.09, is 
the absolute extremum of f(tr) for tr > 0, and f(tr) > 0. Thus, we obtain the exact relation 
between the return time and the initial phase 0 q maximizing Asq, 

sin(t;o + 20*) = l, (16) 

TT 

or equivalently, t* g + 20o = —h 27m, for n > 0. 

Li 
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For n = 1 , this gives <pQ ~ 1.88. Substituting {t*Q,<pQ) into Eq. (14), we finally get the 
maximum change in energy Asq ~ 3.17 Up, where Up = EljAuJ^ or simply 1/4 in rescaled 
units. 

Remarkably, we have found a very simple relation in Eq. (16) between the initial phase 
of the laser and the return time by maximizing the change in energy instead of the return 
energy itself. This relation is borne out by the numerical results for the local maxima of 
x{tr,(p) shown in Table 1 of Ref. [13], when taking care to subtract 7 r /2 from each of the 
phases in that table to reconcile our choice of sin versus cos for the laser field term in 
Hamiltonian ( 2 ). All odd multiples of 7 r /2 appear in that table because local extrema are 
being considered, and /(A) < 0 for the minima, requiring sin(A + 20) = — 1 . Equation (16) 
follows naturally from maximizing As, and in our view, Ae is the fundamental quantity to 
maximize because it is the energy gained by the electron from the laser field. Thus, we will 
only consider maximizing As for the remainder of this paper. 

Now we will consider the effect of allowing Xq 0 on the trajectories that maximize 
As in the SEA, because in the three-step model the electron’s tunnel ionization implies its 
starting position is not exactly xq = 0. Though the consequences of xq 7 ^ 0 are investigated 
in Ref. [13], these simulations include the Coulomb field. Therefore it is not clear whether 
the observed effects are due to the inclusion of the Coulomb field or the relaxation of the 

= 0 condition, and we intend to clarify this. We have already shown that x(t = 0) = 0 
will persist when Xq 0, but to see the effect on the other quantities Ae in Eq. (11) must be 
maximized. If we restrict Xq <^1, then we can treat the effect of Xq 7 ^ 0 on the maximization 
of As perturbatively. This is not an unreasonable restriction because xq just needs to be 
small compared to the quiver radius. For example, if one assumes that ionization occurs 
around 1 a.u., the laser intensity is / = 6 x 10^^ W • cm“^, and the laser wavelength is 800 
nm, the corresponding value of Xq in rescaled units is Xq ~ 0.08. For the same frequency 
but an intensity of / = 10^® W • cm“^, which is used for numerical simulations in this paper, 
the corresponding value of Xq in rescaled units is even smaller: Xq ~ 6 x 10“^. To do the 
calculation, we assume first order perturbation expansions of t* and 0 * in integral powers of 
Xq, 


(17) 

(18) 

6 


t* — t*Q + tlXo + 0(xq), 

0 * = 00 + 012)0 + U^ixl), 





Figure 1. (color online) Dependence on xq of the maximum energy exchange Ae* (upper left panel), 
return time t* (upper right panel), laser phase 0* (lower left panel) and t* + (j)* (lower right panel) 
from Hamiltonian (2). Solid (red) curves are the numerical values and the dashed (blue) lines are 
the first order approximations. Grey dotted lines are the values for xq = 0. xq is in units of E^juP' ^ 
Ae* is in units of Up, t* and (j)* are in radians. 

with the zeroth order terms taking their values from the Xq = 0 case. We carry out the 
calculation to first order by computing {ti, ^i). Recalling Eqs. (13) and (14) and setting the 
partial derivatives of Eq. (11) with respect to tr and (f) equal to 0, we get 


^tr(^^o) + ^0 dtrQ — 0) 

( 19 ) 

^^(Aso) + XQd^g = 0. 

( 20 ) 


Now we substitute the perturbation expansions (17) and (18) into Eqs. (19) and (20) and 
discard terms of order 0{xq). By Taylor-expanding the derivatives of A^o to first order in xq 
and applying phase relation (16) to the (t* 0 ) 0o) terms, we get the system of linear equations 

2A£o\ a a _ (^trg\ 

2A£^ 4Ae*oJ [d,pgj 

where “|o” means to evaluate at the parameters from the xq = 0 case. Equation (21) yields 
(ti,(pi) ^ (0.38,-0.29). 


( 21 ) 
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The results of these approximations are compared against numerical maximization of 
Eq. (11) in Fig. 1. These approximations do work well for |xo| <C 1, and even appear to 
do a fair job closer to |xo| = 1. Note that the shift of the final phase tr + (/> is slightly less 
pronounced than for tr and (j), as its first order coefficient is small, = 0.09, so that 

t* + (p* appears fairly constant. 

We notice on the upper left panel of Fig. 1 that Ae* significantly exceeds the well-known 
value of 3.17 Up. This has already been noted in Ref. [13] [and is also easily seen from 
Eq. (10) using Eq. (9)]: The maximum is expected to be at 8 Up which occurs at = vr 
(e.g., with parameters (p = tt and U = Stt with Pq = 1). However only the region close to the 
core is physically relevant since the electron is initially bound to the core. Therefore we do 
not expect significant variations of initial laser phase, ionization time and maximum return 
energy from the values obtained from xo = 0. In summary, the typical energy-maximizing 
recollision occurs within one laser cycle and brings an amount of energy close to 3 Up. When 
the laser field is large, we expect a similar result by taking into account the Coulomb field, 
and this is what we are going to investigate in the next section. 

When examining the effect of the Coulomb field on the energy-maximizing trajectories, 
we need to study the location of their initial conditions in phase space in detail. Before 
proceeding it is enlightening to look at the phase space picture of energy-maximizing trajec¬ 
tories in the SEA. In order to locate these trajectories in phase space, we consider a Poincare 
section of the trajectories. In the left panel of Fig. 2, we consider a grid of initial conditions 
(xo,Po) beginning on the Poincare section 


X = 0 with X < 0, (22) 

and integrated for six laser cycles, with xq G [— 0 . 08 , 0 .OSjEo/o;^ and po G [— 0 . 2 , 0 . 35 ]Eo/a;. 
We notice that the Poincare section determines the initial phase such that cos (p = —po and 
sin0 > 0 . Since x = pq + cos{t + <p), the Poincare section is a stroboscopic plot with period 
27r in the SEA. The condition i: = 0 is a natural section because we have just shown that 
all of the trajectories maximizing Ae satisfy this condition in the SEA. For all trajectories 
beginning on this Poincare section, the electron begins at rest and starts to move to the 
left. The reversal of the laser field causes the electron to reverse its direction and move to 
the right, towards the core again. It is then possible for a recollision to occur, and thus we 
define a “recollision” as crossing x = 0 moving to the right, i.e. with i: > 0. 




Xo Xo 


Figure 2. (color online) Energy exchange As: at first recollision (after a six laser cycle integration) 
as a function of initial condition, (xo,2^o)- The initial phase 0 is fixed by Poincare section condition 
(22). Darker points correspond to higher energy. Left panel: The initial conditions (xo,po) bringing 
the maximum energy As: at first recollision for each xq are the light grey (blue) curve. Right panel: 
Zoom of the black-framed region on the left. Solid grey (colored) lines separate regions of initial 
conditions that have their first recollision during different laser cycles, with the colors corresponding 
to Fig. 3. Xo is in units of Eo/cj^, po is in units of Eo/cj, As: is in units of Up. 


In Fig. 2, we show the energy As: brought back by each initial condition on its first 
recollision. When starting at the top of the left panel of Fig. 2 and moving down, we observe 
a continuous gradient of As:, until a series of discontinuous changes beginning at small po > 0 
and Xq < 0. These discontinuities are seen clearly in the right panel of Fig. 2, which shows 
that the first discontinuous change in As: is a sudden increase in energy, followed by another 
continuous decreasing gradient of As:. This behavior repeats over and over as po gets closer 
and closer to 0. It turns out that these sudden changes in energy as the initial conditions 
are varied correspond to changes to the laser cycle during which an initial condition has its 
first recollision, referred to as delayed recollisions in what follows as opposed to immediate 
recollisions which occur during the first laser cycle. The delayed recollisions are analogous 
to the “higher-order returns” of Ref. [6]. In Fig. 3 we determine which laser cycle the first 
recollision will occur for a given initial condition analytically using Eq. (4) and the Poincare 
section condition (22). The general structure in Fig. 3 is that immediate recollisions in the 
“1” region are separated from the immediate ionizations in the “0” region by a small triangle 
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Figure 3. (color online) Number of iterations of Poincare map V until the first recollision has 
occurred, as a function of initial conditions (a:o,Po)) where cj) is fixed by the Poincare section 
condition (22). In the SFA, a single iteration of this map always corresponds exactly to the 


duration of one laser cycle, xq is in units of po is in units of EqIuj. 


in which the delayed recollisions live. A closer look at the delayed recollision triangle in the 


inset of Fig. 3 shows that the triangle is bounded below by the line po = 0 and its rightmost 


point is (xo,po) = (0,0). Furthermore, the triangle is clearly stratified, with each stratum 
corresponding to a delayed recollision with a particular delay. The strata are well ordered in 
phase space: as po decreases, the delay time increases. These observations are explained by 
considering the mechanical scenario in this region: All trajectories starting on the section 
will first move to the left, and as is evident from Eq. (4), they will eventually come to 
rest (at approximately half of the laser cycle for small po), move to the right, and come to 
rest again, this time on the Poincare section (with exactly one laser cycle elapsed). They 
will have a displacement from their initial position of Ax = 2 'kpq. Therefore, trajectories 
beginning with Xq < 0 and a small enough po > 0 such that Xq + Ax < 0 will not experience 
a recollision in the first laser cycle: they must wait for enough laser cycles n such that the 
displacements accumulate and satisfy xq + nAx > 0. The smaller the po, the more laser 
cycles are needed for a recollision to occur, which explains the ordering of the strata. This 
argument also explains why there are no recollisions if xq < 0 and po < 0, and none for xq > 0 
and Po < — Xo/27r. The lines that form the boundaries between the delayed recollisions of 
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differing delay times are plotted in the right panel of Fig. 2. We note that the boundaries 
occur precisely where the return energy has discontinuous jumps, as claimed earlier. 

It should be noted that the highest energies associated with delayed recollisions are an 
order of magnitude smaller than the immediate recollisions, which is reflected in the change 
of color scale between the left and right panels of Fig. 2. Thus, when searching for energy- 
maximizing trajectories in the SFA, these kinds of trajectories appear irrelevant, not bringing 
enough energy to play a significant role in recollision-driven processes. However it should 
be emphasized that this is only on their first recollision. In the SFA, these trajectories with 
very small po will have many recollisions, as is evident from Eq. (4). With x = po+cos{t+(f)), 
the return kinetic energy of the later recollisions can only get slightly higher than 2 Up, still 
significantly below the usual prediction of 3.17 Up. Moreover, the maximum return kinetic 
energy of any recollision (i.e. not necessarily the trajectory’s first recollision) occurring after 
the first laser cycle for a trajectory initiated near a; = 0 is bounded at approximately 
2.4 Up [13]. This changes when we turn on the Coulomb field: even if the laser field is large, 
the Couloumb field allows delayed recollisions to bring back energy greater than 3.17 Up on 
their very first recollision. In addition, the trajectories experiencing delayed recollision are 
particularly interesting because they spend a longer time close to the core where they can 
exchange energy and experience the influence of the core potential. 


III. MAXIMIZING THE RETURN ENERGY WITH THE COULOMB FIELD: 
NONLINEAR DYNAMICS OF DELAYED RECOLLISIONS 

A. Delayed Recollisions 

In Sec. II, we have identified the effects of Xq 7 ^ 0 on the recollision mechanism in the 
strong field approximation. In particular we have assessed the effect of Xq 7 ^ 0 on the maxi¬ 
mum Ae trajectories. We have shown that the energy-maximizing trajectories are the ones 
which recollide within one laser cycle. In this section, we turn on the Coulomb interac¬ 
tion and investigate how the strong field approximation picture for recollision trajectories is 
affected by this interaction. The Hamiltonian, in atomic units, now becomes 

H {x, p, t) = I-(p+—cos{ut+ (/))] I , (23) 

z\uj / +1 
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where we have used a soft-Coulomb potential with a softening parameter of 1 [4, 14], The 
general phase space picture associated with Hamiltonian (23) is composed of two distinct 
regions [15]: the region near the core, referred to as the bounded region, where the Coulomb 
attraction dominates compared to the laser forcing, and the region sufficiently far from the 
core, referred to as the unbounded region, where the laser held dominates. For low enough 
laser intensities, the electron never escapes the bounded region. A common assumption in the 
unbounded region is to neglect the Coulomb interaction for sufficiently large intensity, which 
corresponds to the second step in the three-step model of recollisions. Here we investigate 
the validity of this assumption, using numerical simulations of the equations of motion 
arising from Hamiltonian (23). In order to focus on the effect of the Coulomb held in 
the unbounded region, we choose the parameters in such a way that the bounded region 
is completely suppressed. In particular, we use a laser intensity I = 10^® W • cm“^ and 
frequency u ~ 0.057 a.u. (corresponding to an 800 nm wavelength). At this laser frequency, 
the bounded region only exists for intensities below / = 4.9 x 10^^ W • cm^. Thus, our choice 
of laser intensity is high enough that almost every initial condition will lead eventually to 
ionization, meaning there are almost no trajectories that remain indehnitely bounded to the 
core. In fact, this intensity is so high that there is not even a potential barrier when the 
laser held is at its maximum, as shown in Fig. 4, so classical ionization can account for all 
ionizations; it is not necessary to invoke a tunneling argument. 

We rescale this Hamiltonian as in Sec. H, yielding 

'H(x,p,t) = J(p + cos(t + (/)))- , / (24) 

where a = /Eq 6.1 x 10“^ and e = /Eq ~ 6.9 x 10“^. The parameter e is the 
effective strength of the Coulomb potential in the presence of the laser held: as the intensity 
goes to inhnity, e goes to zero as does e/a so the Coulomb potential becomes negligible. In 
this parameter regime it would appear that the Coulomb potential is ignorable due to the 
small magnitude of the maximum Coulomb energy, e/a ~ 0.046 Up. However, the relative 
importance of terms in the Hamiltonian is not determined by the absolute values of the 
terms themselves but by the absolute values of the gradients of the terms, because it is the 
gradients that actually appear in the equations of motion. In this parameter regime, the 
Coulomb force, equal to —ex/{x‘^ + a^Y^‘^ has a maximum amplitude of 2e/(3a^\/3) ~ 0.7Eq, 
which is comparable to the maximum electric force, and therefore cannot be neglected. 
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Figure 4. The potential of Hamiltonian (24) in the length gauge, i.e. V(x) + xE{t), as a function of 
X. Grey dotted is with E{t) = 0, black dashed with E{t) = 0.25 £"0 (the maximum field amplitude 
when I = 6.3 x 10^^ W • cm“^), and solid black with E{t) = Eq. x is in units Eojoj'^ and V is in 
units of Up. 


In the same way as we did in Sec. II, we examine the trajectories that maximize the 
change in energy Ae between some starting position at t = 0 and the final position at 
the center of the core, a; = 0, at the return time tr- With the inclusion of the Coulomb field, 
the translational invariance of the SFA is broken, meaning the momentum p is no longer 
conserved. Thus, this system is no longer integrable, so trajectories cannot be obtained 
analytically. Here we numerically integrate the equations of motion associated with Hamil¬ 
tonian (24) in order to find the trajectories that maximize Ae. We notice that the discrete 
symmetry (3) remains, meaning we may still confine our simulations to initial conditions 
with (f) G [0,7r). Using a symplectic integrator [16], we again consider a large number of 
initial conditions (typically of the order of one thousand for each value of xq) beginning on 
the Poincare section (22), and look for the trajectories that yield the highest As at the first 
recollision. The range of initial conditions is the same as in the previous section, and with 
this choice of laser parameters the ranges of initial conditions in atomic units correspond to 
Xq ~ [—13,13] a.u. and po ~ [—1.9,3.3] a.u. 

In the top left panel of Fig. 5, we show the energy As brought back by each initial 
condition on its first recollision. Comparing this with Fig. 2, we can evaluate the effect of 
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Figure 5. (color online) Energy exchange A^: at first recollision after a six laser cycle integration of 

Hamiltonian (24) with a ^ Q.lx 10“^ and e ^ 6.9 x 10“^, as a function of initial conditions, (xq^pq) 

with the initial phase (f) is fixed by Poincare section condition (22). Darker points correspond to 

higher energy. Top left panel: The initial conditions (xo,po) bringing the maximum energy As: at 

first recollision in the first laser cycle for each xq are the light grey (blue) curve. Top right panel: 

Same as top left, with the stable manifold of O plotted in black and its unstable manifold 

in grey. Bottom left panel: Zoom of the lower grey (magenta) framed region on top. Grey (colored) 

curves are the initial conditions (xo,po) bringing the maximum energy at first recollision for each 

Xq and iteration of Poincare map V during which the recollision happens. The colors correspond 

to Figs. 3 and 6. Bottom right panel: Zoom of upper grey (green) framed region above, xq is in 
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units of E^jbj^^ po is in units of E^juj^ As: is in units of Up. 

































the Coulomb field on the recoil!ding trajectories. The shape of the curve of initial conditions 
leading to Ae-maximizing immediate recollisions is deformed compared to the SFA, due to 
the presence of the Coulomb field. The actual energies of these recollisions as a function of Xq 
are seen more clearly in Fig. 7. We see that the immediate recollision curve is mostly close to 
the SFA curve, especially as |xo| increases. At such a high laser intensity, its deviations from 
the SFA curve are in large part due to the leading Coulomb correction to the final return 
energy e(tr) [17], which is indeed even closer to the SFA curve. Therefore, we now see clearly 
the differing effects of varying Xq versus turning on the Coulomb field on the return energy. 
If one considers the total energy at return (as opposed to the energy difference or only the 
return kinetic energy), then the effect of the Coulomb field on the value of the return energy 
compared with the SFA is small for this laser intensity. The major deviations from 3.17 Up 
come from varying xq, and in the previous section we account for these effects using the 
SFA. 

Now we reexamine the delayed recollisions with the Coulomb field on. The continuous 
gradient of Ae in the upper left panel of Fig. 5 near the maximum energy immediate recolli¬ 
sions is again interrupted, at low values of po) by a series of jumps, seen more closely in the 
bottom left panel of Fig. 5. This is where the delayed recollisions originate, which can be 
confirmed by inspecting Fig. 6. Like Fig. 3, Fig. 6 shows the iteration of the Poincare map 
during which each initial condition has its first recollision. The Poincare map is the discrete 
map (a;j+i,pj+i) = V{xi,pi) that takes a point {xi,pi) on the Poincare section (22) to the 
next point (xi+i,Pi^i) at which the trajectory started with {xi,Pi) pierces the surface of sec¬ 
tion. Thus, for the range of initial conditions we consider, the iteration of the map V during 
which an initial condition {xq,pq) has its first recollision is the smallest n such that Xn > 0. 
As we stated earlier, in the SFA the map V corresponds exactly to the duration of one laser 
cycle. With the Coulomb field on and a large laser intensity, this correspondence continues 
to hold only approximately. Still, for the remainder of the paper, when we refer to a “cycle” 
during which a recollision occurs, we really mean the iteration n of the Poincare map V. 
Comparing Fig. 6 and Fig. 3 makes the significant qualitative differences between the SFA 
and the full Hamiltonian (24) apparent. For example, the Coulomb focusing effect is clearly 
manifested by the greater area of phase space that leads to recollisions compared with the 
SFA picture. The curve that separates immediate ionizations from recolliding trajectories 
in Fig. 3 is plotted on the upper panel of Fig. 6. We see that a substantial area of the im- 
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mediate ionizations in the SFA are converted into recolliding trajectories with the Coulomb 
held on, especially in the core region (small Xq). In particular, the delayed recollisions are 
enhanced by the Coulomb held. In the previous section we observe from Fig. 3 that the 
delayed recollision initial conditions lay between the immediate recollisions and immediate 
ionizations in the region of phase space with a small po > 0 and xq < 0. With the Coulomb 
held on, we see in Fig. 6 that the delayed recollision boundary extends to xq > 0 and near 
Xo = 0 it takes a signihcant dip in the direction po < 0, so the presence of the Coulomb held 
enables more delayed recollisions. 

In addition, the Coulomb held allows certain delayed recollisions to bring energy that 
slightly exceeds the well known 3.17 Up upon their very hrst recollision. This is shown in 
the lower left panel of Fig. 5, where although we have changed the distribution of colors on 
the color scale compared with the panel above it, the overall scale remains the same. It is 
true that most of each recollision region of a particular delay contains trajectories that hrst 
recollide with a low amount of energy close to that brought by delayed recollisions in the SFA, 
^ 0.5 Up. However, the energy at hrst recollision rapidly increases near the boundaries of 
the recollision regions of differing delays. The initial conditions bringing back the maximum 
energy for each Xq and delay time are the colored curves. Actually these curves are made up 
of discrete points (one for each xq on our grid) that in most places are so close together that 
they appear to form a continuous curve. The energies of the second-cycle delayed recollisions 
as a function of xq are plotted in Fig. 7. We observe that for any xq, there exists a second- 
cycle delayed recollision that brings back energy exceeding 3.17 Up. Moreover, the amount 
of energy brought back by these delayed recollisions does not display the approximately 
linear dependence of Ae^ax on xq characteristic of the immediate recollisions, but rather 
approaches a constant value for large |xo| and has a jump near x = 0. Surprisingly, the 
Coulomb well is all but invisible at maximum laser amplitude, yet its presence causes certain 
delayed recollisions to recollide with an energy above the usual SFA prediction on their very 
first recollision. Though we do not show the maximum energies brought by higher order 
delayed recollisions in Fig. 7, we argue in the next section that for every Xq in the range we 
consider, there is a recollision of arbitrary delay that carries the same energy as we observe 
the maximum energy second-cycle recollisions carry [18]. Therefore, during every laser cycle 
there are energetic recollisions with energies close to 3.17 Up, and this is due to the Coulomb 
field. We emphasize that the SFA misses this completely - as discussed in Section II, the 
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Figure 6. (color online) Number of iterations of Poincare map V until the first recollision has 
occured after a six laser cycle integration, as a function of initial condition (xo,po)- The initial 
phase (j) is fixed by Poincare section condition (22). Top left panel: The black dotted line separates 
immediately ionizing trajectories from recolliding trajectories in the SFA. Top right panel: Same as 
top left panel, with stable manifold of O plotted in black. Location of O- marked as the grey 
(red) cross. Bottom: Zoom of the grey (magenta) framed region above, xq is in units of ^ po 

is in units of 


energetic immediate recollisions are ionized after one laser cycle, and only trajectories with 
sufficiently low drift momenta will continue to recollide with energies ^ 2.4 In the next 
section, we explain the nonlinear dynamical origin of the delayed recollisions, and the origin 


17 

























Figure 7. (color online) Maximum change in energy As as a function of xq. Thick black line: SFA, 
immediate recollision. Grey (light blue) curve: Coulomb included, immediate recollision. Light 
grey (cyan points): Coulomb included, second cycle delayed recollision. Thin black (blue) curve: 
Coulomb included, immediate recollision final energy, s(t^). Dotted black line: 3.17 Up. Dotted 
grey (red) line: Estimate of delayed recollision energy “cutoff.” xq is in units of E^juj^ e is in 
units of Up. 

of the cutoff in their energy exchange. 

B. Building Blocks of Recollisions: Periodic Orbits and Their Invariant Manifolds 

The qualitative character of both immediate and delayed recolliding trajectories can be 
understood by considering some periodic orbits of the system. The most relevant periodic 
orbits are O and (9_, shown as the insets on the left plot of Fig. 8. The importance of O 
for recollisions was first described in Ref. [17], and in Ref. [19] it is shown that (9_ (and 
(9+, its symmetric copy via symmetry (3)) underlie ionization stabilization in ultra-intense 
laser fields. Due to our choice of Poincare section and our interest in initial conditions 
beginning near the core, we need not consider here. Each orbit has a period of exactly 
one laser cycle. For the parameters we have chosen, both O and O- are hyperbolic, i.e. they 
are unstable. Thus, initial conditions in the vicinity of these orbits may take trajectories 
qualitatively similar to the orbits for some time, before eventually ionizing. 

This is precisely what we observe when we examine both the immediate and delayed 
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Figure 8. (color online) Left: Two typical Ae-maximizing trajectories in x-p phase space. An 
immediate recollision is in black (blue), and a third-cycle delayed recollision is in grey (red). Left 
inset: Periodic orbit O. Right inset: Periodic orbit 0-. Right: Ae-maximizing trajectories of the 
left figure plotted as x vs. t, upper, and p vs. t, lower. The immediate recollision is also shown 
as a dotted black (blue) line, translated in time so its laser phase matches that of the delayed 
recollision. x is in units of Eojuj'^, p is in units of Eq/lo, t is in radians. 


recollision trajectories in phase space and compare them with the periodic orbits. In the left 
panel of Fig. 8, we have plotted an energetic immediate recollision and an energetic third- 
cycle delayed recollision. The immediate recollision resembles the motion of O, as was noted 
in Ref.[17]. O's main characteristics are motion between x ~ —Eq/uj'^ and x ~ with 

a small momentum kick due to the Coulomb field upon recollision. Correspondingly, the 
energetic immediate recollision also travels out to a: ~ —Eq/uj"^, before recolliding, receiving 
a small momentum kick, and ionizing forever. Delayed recollisions on the other hand spend 
some time following the orbit O -, before reaching the part of phase space where immediate 
recollisions live and the motion of O dominates. O- is distinguished by its motion between 
X ~ —2Eq/u}‘^ and a: ~ 0, with a comparatively large momentum kick at the core. Inspection 
of the third-cycle delayed recollision confirms that it has exactly this behavior for two laser 
cycles. After the end of the second laser cycle, the momentum kick received by the electron is 
so large that the trajectory is moved into the region of phase space from which the energetic 
immediate recollisions originate. Thus, delayed recollisions follow the motion of O- for some 


19 

































number of laser cycles until they are converted into immediate recollisions. At this stage, 
their motion is the same as an immediate recollision. This can be observed from the right 
panels of Fig. 8, where we have plotted the energetic recollisions of the left panel as both x 
vs. t and p vs. t, with the immediate recollision also translated forward in time so that its 
laser phase matches that of the delayed recollision. Clearly, these trajectories nearly overlap, 
illustrating that the delayed recollisions eventually become immediate recollisions. 

The arrangement of the delayed recollision initial conditions in phase space and their 
transport to the realm of immediate recollisions is determined by the stable and unstable 
manifolds, respectively, of O. The stable manifold W® of a periodic orbit O is defined as 
the set of initial conditions that approach (9 as t —?• oo. Conversely, the unstable manifold 
>V“ is the set of initial conditions that approach O as t ^ —oo, meaning those trajectories 
escape from the orbit moving forward in time. It was already shown in Ref. [17] that O’s 
invariant manifolds regulate the recollision dynamics. There, the initial conditions leading 
to many recollisions were found to be concentrated near W®, and trajectories would follow 
>V“ on their way to ionization. It turns out that these manifolds also organize the delayed 
recollisions: delayed recollisions begin near O’s stable manifold and then follow its unstable 
manifold to eventually recollide with the core. The former is seen in the right panel of 
Fig. 6, where we show W® [20]. TV® is obtained from >V“, which we computed using the 
algorithm of Ref. [21], by the time-reversal symmetry {x,p, (p) {x, —p, tt — 0). Comparing 

the left and right upper panels of Fig. 6 shows that all the delayed recollisions come from a 
region of phase space in which the curves of the stable manifold are dense (compared with 
other regions of phase space). Looking at the delayed recollision region more closely in the 
lower panel of Fig. 6 shows that the stable manifold actually separates regions with differing 
recollision delay times. This is seen at the boundaries between immediate/second cycle 
recollisions, second/third cycle recollisions, and third/fourth cycle recollisions. We expect 
that if the manifold is computed for longer times then it would also be seen to separate 
the higher order recollisions. We note that there are parts of the manifold that are entirely 
contained in a region of a particular recollision delay time. While these parts of the manifold 
are not relevant for separating initial conditions with different delay times, we expect they 
separate trajectories with another type of qualitative difference; for example, they may 
separate regions of initial conditions that have differing total numbers of recollisions, as in 
Ref. [17]. Additionally, the delayed recollisions are indeed located near the orbit 0-, as 


20 




X 

Figure 9. (color online) V{xo^po) for the initial conditions (xq^pq) in Fig.6 leading to a second-cycle 
delayed recollision. Color scale shows xq. Thick grey (light blue) curve is the initial conditions 
leading to a maximum As: immediate recollision. Thin grey curve is of O . Inset: Zoom of the 
grey (magenta)-framed region, x is in units of Eo/(jj‘^^ p is in units of Eq/lj. 

was suggested earlier by noting the similarities of the delayed recollision trajectories to the 
periodic orbit. 

When the initial conditions leading to delayed recollision are integrated, they follow 
the unstable manifold of (9, and this process underlies our observations regarding the maxi¬ 
mum energy available to a delayed recollision and the particular dependence of the maximum 
energy on xq. Consider an initial condition leading to an n-cycle delayed recollision. By 
definition, after one iteration of the Poincare map, it must move into the region of phase 
space from which (n — l)-cycle recollisions originate, with x < 0. This will happen repeat¬ 
edly until the iteration immediately prior to recollision, when the trajectory reaches the 
immediate recollision region of phase space. Additionally, the trajectories accomplish this 
motion by following VF^. Therefore, the recollision energies accessible to delayed recollisions 
are the immediate recollision energies in the vicinity of the unstable manifold. These can 
be read off by looking at the upper right panel of Fig. 5. VF^ enters the immediate recol¬ 
lision region above the region where VF^ is dense and delayed recollisions proliferate. We 
observe that in this region VF^ is thickest for smaller po? in this region the recollision 
energies are low, ^ 0.5 Up. This explains why most of the delayed recollisions recollide with 
an energy ^ 0.5 t/p, as we remarked earlier in reference to the lower left panel of Fig. 5. 
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However, W“ does enter the region of high energy immediate recollisions, and in fact it is 
seen to intersect the curve of maximum Ae immediate recollisions in the lower right panel 
of Fig. 5. The delayed recollisions that end up in the region near this intersection are the 
delayed recollisions maximizing Ae. 

We look more carefully at the intersection of >V“ with the curve of Ae-maximizing imme¬ 
diate recollisions in Fig. 9. Here we have also plotted many second-cycle recollisions on their 
hrst return to Poincare section (22). Firstly, the points are seen clearly to lie very close to 
>V“, confirming our earlier claim that the delayed recollisions follow W“. Next, we observe 
that there are always continuous gradients of Xq perpendicular to the manifold, including 
in the region where the manifold intersects the curve of maximum energy immediate recol¬ 
lisions. This region is magnified in the inset of Fig. 9. This is exactly why for any xq, it is 
possible to find a second-cycle delayed recollision with As exceeding 3.17 Up. The energy of 
these delayed recollisions is 

ASmax A£ijiax,imni(^p) T AV(^Xp^ XqJ ^ (2^) 

where A£max,imm(3:) is the maximum energy of an immediate recollision beginning at x, Xp 
is the position of the delayed recollision’s return to the Poincare section on the cycle imme¬ 
diately prior to its recollision, and AH {xp, xq) is the change in energy between the delayed 
recollision at t = 0 and the time it returned to the Poincare section prior to recollision. 
Because the section condition is zero velocity, there is no change in kinetic energy so AH is 
just the change in Coulomb potential energy. Therefore we can obtain an estimate of the 
minimal A^max for a delayed recollision by taking AH(xp,xo) = —eja and A£max,imm to be 
the minimal maximum energy accessible to an immediate recollision in the vicinity of W“. 
This estimate is plotted for second-cycle delayed recollisions in Fig. 7, and does a good job of 
providing a “cutoff” for the minimal maximum energy second-cycle delayed recollision. The 
remaining shape of the maximum energies for the delayed recollisions resembles the Coulomb 
potential and is mostly due to AH {xp, Xq). Intuitively, because Xp takes a quite limited range 
of values, A£:inax,imm(a^p) is approximately constant and AH(xp,xo) ~ V{xp) — H(xo), where 
Xp is some position where TV" intersects the maximum energy immediate recollision curve. 
This explains why the delayed recollision energy curve approximately takes the shape of 
minus the Coulomb potential. 

In Fig. 9, is seen to intersect the curve of Ae-maximizing immediate recollisions in two 
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Figure 10. (color online) Energy exchange Ae at first recollision as a function of initial conditions 
{xq,pq). a region of second-cycle delayed recollisions is highly magnified here, revealing two local 
maxima in Ae. Second-cycle delayed recollision initial conditions are in the colored region, imme¬ 
diate recollisions in the white region. Darker colors correspond to higher recollision energies, xq is 
in units of EqIuP', pq is in units of Eq/w, Ae is in units of Up. 

places. In between the two intersections, crosses into a region of lower energy recollisions, 
seen in the lower right panel of Fig. 5, and thus we should expect that there should actually 
be two energy maximizing delayed recollisions for each Xq. Indeed this is what we see, and 
this is the cause for the thickness of the set of points giving energy maximizing delayed 
recollisions in Fig. 7. This is seen even more clearly in Fig. 10, where we have magnified a 
region containing second-cycle delayed recollision initial conditions and shown their energy 
exchange upon their first recollisions. We observe two ridges, corresponding to the predicted 
two maximum energy recollisions for each Xq. 

It is also possible to find higher order delayed recollisions of an arbitrary delay that return 
with an energy at least at the delayed recollision cutoff. As we have argued, every laser cycle 
a delayed n-cycle recollision will move into the region of phase space from which the n — 1- 
cycle recollisions originate, with a: < 0. Therefore, eventually any n-cycle delayed recollision 
will end up in the second-cycle delayed recollision region of phase space with a: < 0, i.e. 
two laser cycles before its recollision. But we have just shown that for any xq, there are 
second cycle delayed recollisions that recollide with an energy above the delayed recollision 
cutoff. Thus, every laser cycle, there are delayed recollisions arriving that bring an energy 
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near 3.17 Up to the core. Though the SFA does a reasonable job of providing the energy 
cutoff, these kinds of trajectories are absent in the SFA, and are a direct consequence of the 
ion’s Coulomb field and the organization of the dynamics by the periodic orbits O and 0± 
and the invariant manifolds of O. 


CONCLUSION 

In summary, we have shown that even in the high laser intensity regime where the strong 
field approximation is expected to hold, the Coulomb field significantly impacts the dynamics 
by allowing delayed recollisions to bring high energy to the core region. We reported the 
effectiveness of Coulomb focusing in a one-dimensional model, which is usually attributed 
to higher dimensional models. We found that while the SFA gives adequate estimates of 
the maximum possible return energy in the strong field regime, it misses the behavior of the 
electron for times beyond the first laser cycle, when the Coulomb field causes trajectories 
to continue to recollide with energies near 3.17 Up. We have unraveled the dynamical 
organization of these delayed recollisions by looking at specific periodic orbits and their 
invariant manifolds. The delayed recollision trajectories that we focused on are important 
not only because of the high energy they may recollide with, but also because they spend 
potentially many laser cycles near the core. Thus, they have many opportunities to exchange 
energy with the ion and the electromagnetic field (due to the acceleration imparted on such 
electrons by the Coulomb force). We expect that the delayed recollisions are the mechanism 
underlying HHG and thus explain the observation of the appearance of the plateau and 
~ 3.17 Up high-harmonic cutoff only after the first laser cycle, as in Ref. [17]. 
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